Improved Kl — > 7r ± e T z/ Form Factor and Phase Space Integral 
with Reduced Model Uncertainty 
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Using the published KTeV sample of 2 million Kl — > ix^e^v decays and a new form factor 
expansion with a rigorous bound on higher order terms @, we present a new determination of the 
Kl — » ir^e^v form factor and phase space integral. Compared to the previous KTeV result 0, the 
uncertainty in the new form factor expansion is negligible and results in an overall uncertainty in 
the phase space integral (Ik) that is a factor of two smaller: I C K — 0.15392 ± 0.00048. 

PACS numbers: 13.25.Es, 14.40.Aq 



Our previous analysis of semileptonic form factors 
contributed to improving the precision in the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix element \V US \ 0,0, 
5] . The published form factor results were based on lin- 
ear, quadratic and pole models: 
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where t = (Pk — Pit) 2 is the expansion parameter as a 
function of the kaon and pion four- momenta, f+(t) = 
/+(*)//+ (0) is the normalized form factor, and A+, A' 
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My are the form factor parameters measured in 
to data. Figure Q shows the i-distribution for 
* ir^e^v phase space with and without a form fac- 
tor. The best fit to the KTeV data, with x 2 /dof = 62/65, 
@ was obtained with fl uad (t) (Eq.EJ. The pole model, 
with one free parameter, also gives a good fit with 
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FIG. 1: calculated t-distribution, dT/dt, for Kl — > ^e^v 
phase space (dashed), and for phase space modified by the 
form factor in (solid line). 



X 2 /dof = 66/66. The linear model, with x 2 /dof — 
81/66, is disfavored. Although the data are well de- 
scribed by the quadratic and pole models, the corre- 
sponding phase space integrals (7^) differ by 0.7%, which 
is much larger than the experimental precision of 0.3%. 
This 0.7% difference was arbitrarily included in the 1^ 
uncertainty for the determination of \V US \. 

The form factor models in Eqs. Ill.il are theoretically 
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flawed: f+° le (t) makes the unlikely assumption that there 
is only one pole, and truncated Taylor expansions vi- 
olate unitarity. These issues are addressed in a new 
parametrization based on analyticity and bounds from 
t — > Kirv decays |2(. In this paper, we present an up- 
dated form factor and phase space integral based on our 
published data and this new parametrization. 
The new parametrization expands in the variable 

_ *fb+ — t — — to 



z(t,t ) 



Vt^t + yjt+ - t 



(4) 



where t + = [Mk + M-n) 2 and to is an arbitrary an- 
chor point. Note that our previous expansion parameter, 
t/M%, has a range of values from zero to about six; the 
new expansion parameter is limited to \z\ < 0.2. The 
form factor expansion in this new variable is 

FUt) = FUto)%^^Y / a k (to,Q 2 )z(t,t ) k , (5) 



<f>(t,h,Q 2 ) 



fe=0 



where do = 1, ate (k > 1) are coefficients to be deter- 
mined, Q 2 is the invariant momentum-squared of the cur- 
rent in the operator product expansion, and the function 
d> is defined as 
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with t± = (Mk±Mk) 2 ■ The advantage of this expansion 
is that higher order terms are bounded by 



A 2 



k=0 
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Using the previously published KTeV sample of 2 mil- 
lion Kl — > 7r ± e =F i/ decays, the smallest overall error on 
Tjf is obtained from a second order "quadratic-z" fit to 
F+(t) in Eq.0 The expansion coefficients ai(t ,Q 2 ) and 
02(^0, Q 2 ) depend on the choice of to and Q 2 ; however, 
the function (t) , the calculated 1^ value, and the un- 
certainty on Iff are independent of to , Q 2 Q ■ We choose 
to fit with Q 2 = 2 GeV 2 /c 2 and t = 0.39i_ because this 
results in no correlation between a\ and a 2 if the experi- 
mental acceptance of Kl — > 7r ± e =F ^ decays is uniform as 
a function of t. 

We fit the KTeV Kl — > n^e^u sample to the normal- 
ized form factor, 



mt) = Fl{t)/FUQ) 



(8) 



The results are 



a 1 


= 1.023 ±0.028(stat)±0.029(syst) 


(9) 




= 1.023 ±0.040 


(10) 


«2 


= 0.75±1.58(stat)±lA7(syst) 


(11) 




= 0.75 ±2.16 


(12) 


X 2 /dof 


= 62.3/65 


(13) 


P12 


= -0.064 . 


(14) 



The second-order term 02 is consistent with zero, and 
the small correlation (pi 2 = —0.064) is due to the 
slightly non-uniform acceptance of the KTeV detector 
0| . The systematic uncertainties are obtained by scaling 
the ratio of systematic-to-statistical uncertainties for the 
quadratic-i model in Table I of pj, 

1 \ 1 \ . . cr syst(A' + , A" ) 

a S y S t(a 1 ,a 2 ) = <x stat (ai, a 2 ) x —7 — -77- . (15) 

Cstat(A_|_, A + j 

Setting at = for k > 3, the resulting phase space 
integral is 



r 



h - 0.15392 ±0.00035(stat)±0.00033(sysi) 



± 0.00006(</i) 
0.15392 ±0.00048 , 



(16) 
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where the uncertainties are statistical, experimental sys- 
tematic, and theoretical modeling of the form factor. The 
theory error of ±0.00006(t/i) is estimated by substituting 
the maximum allowed value for 03 as follows, 



03 = »J A 2 ^ — clq 2 — min(ai 2 ) — min(a 2 2 ) (18) 
~ V170 - 1 - 1 - 0- 13 , (19) 



and then repeating the fit for a\ and a 2 . Note that a lin- 
ear fit in z gives a statistical error on 1% that is almost a 
factor of two smaller than the error from the second-order 
fit 0- However, the linear fit relies on Sa 2 ~ 13 from the 
bound in Eq. and the corresponding uncertainty in 7^- 
is ~ 0.0023, more than four times larger than the overall 
error from the second order fit (Eq. 1171) . 

Comparing to our previous results in 1], our new 
central value for 1^ lies between the 1^ values de- 
termined from the quadratic^ and pole models. The 
new experimental error (statistical plus systematic), 
S(Ik) = 0.00048, is slightly larger than the experimen- 
tal error reported in [| because both Kl — > n ± e T i' and 
Kl—* were used to determine the form factors 

in 0], while only the K e3 mode is considered in this up- 
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dated result. The overall error of 5I e K = 0.00048 (Ea.H7jl 
is two times smaller than our previously reported error; 
this reduced error is achieved by reducing the form fac- 
tor model uncertainty from ±0.00095 in [l| to a negligible 
0.00006. 

Figure [3 shows the dT/dt distribution for different 
models used to fit the KTeV data. In all cases the dis- 
tribution is divided by (dT/dt) from the linear model 
(Eq. with A+ = 0.0283) so that deviations from the 
linear model are more clearly visible. Our quadratic- z fit 
(Eq. [SJ is much closer to our previous best result with a 
quadratic-i model (Eq. |2J) . A simple parametrization of 
the i-dependence for our fit quadratic- z model is 



(rfr/dt) quad ~ z 

(0.0020286)* 2 
where i = t/M 2 . 



(dT/dt) lin x [1 - (0.01094)4- 
(0.78398 x 10" 4 )i 3 ] , 



(20) 
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FIG. 2: Calculated phase space ratios versus t/M 2 , 
{dT/dt) model /{dT/dt) l "\ where model = pole (Eq. 0, 
quadratic- 1 (Eq. |5J and quadratic- z (Eq. 0. The reference 
model (dT/dt) lin uses f l j n (t) = 1 + 0.0283t/M*. 
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